Signatures of topological phase transitions in mesoscopic superconducting rings 
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We investigate Josephson currents in mesoscopic rings with a weak link which are in or near a 
topological superconducting phase. As a paradigmatic example, we consider the Kitaev model of a 
spinless p-wave superconductor in one dimension, emphasizing how this model emerges from more 
realistic settings based on semiconductor nanowires. We show that the flux periodicity of the Joseph- 
son current provides signatures of the topological phase transition and the emergence of Majorana 
fermions situated on both sides of the weak link even when fermion parity is not a good quantum 
number. In large rings, the Majorana fermions hybridize only across the weak link. In this case, 
the Josephson current is h/e periodic in the flux threading the loop when fermion parity is a good 
quantum number but reverts to the more conventional h/2e periodicity in the presence of fermion- 
parity changing relaxation processes. In mesoscopic rings, the Majorana fermions also hybridize 
through their overlap in the interior of the superconducting ring. We find that in the topological 
superconducting phase, this gives rise to an /i/e-periodic contribution even when fermion parity is 
not conserved and that this contribution exhibits a peak near the topological phase transition. This 
signature of the topological phase transition is robust to the effects of disorder. As a byproduct, 
we find that close to the topological phase transition, disorder drives the system deeper into the 
topological phase. This is in stark contrast to the known behavior far from the phase transition, 
where disorder tends to suppress the topological phase. 



I. INTRODUCTION 

On the way to scalable quantum information process- 
ing Majorana fermions (MF) in topological supercon- 
ductors are a promising candidate for the implemen- 
tation of quantum bits in solid-state devices Since 
information in such systems is stored and processed 
in a nonlocal fashion by means of their non-Abelian 
statistics Major ana-based qubits are immune to local 
fermionic parity conserving perturbations which impair 
other qubit realizations. Manipulation of such topolog- 
ically protected qubits requires braiding of MFs, which 
is well-defined only in two dimensions. However, differ- 
ent schemes have been proposecP^ to enable braiding of 
MF in one dimensional systems. Several suggestions for 
one-dimensional physical realizations that host Majorana 
bound states (MBS) have been made. These suggestions 
are based on conventional superconductors in proxim- 
ity to various systems including a topological insulator 
edge,^ sem icond uctor wires in a magnetic field and 
half metals!™ 

Recently, more realistic investigations have elucidat- 
ed the effects of interactions and disorder. In gen- 
eral interactions^^ or disorder with short range cor- 
relations^^ can greatly affect the range of parameters 
in which the system supports topological boundary states 
or even cause the topological phase to break down com- 
pletely if the interaction^ or disorder^ strength exceed 
certain critical values. Long-range correlated disorder in 
a topological superconductor creates nontopological do- 
mains with MF localized at the domain walls pSHSU g ev _ 
eral proposals have been put forward to ac cess MFs ex- 
perimentally based on interferometr}PIl2Zl or tr ansport 
properties such as tunneling conductance peakp^SOU 



half- intege r con ductance plateaus,^ or signatures in the 
shot noise! 31 ! 32 ! Recent experiments have reported pos- 
sible signatures of Majorana bound states in the differ- 
ential tunneling conductance of semiconductor quantum 
w j res j33HHU 

A more specific way of detecting MFs is to measure 
the Josephson current acro ss a weak link between two 
topological superconductors^ 6 1 9 * 10 1 36 1 that arises due to 
a phase difference of the superconducting order param- 
eters. If the weak link is incorporated into a ring made 
of a conventional superconductor, the current flowing 
through the junction is a periodic function of flux with 
period h/2e (corresponding to 2n periodicity in the phase 
difference), associated with the transfer of Cooper pairs 
across the junction. In a ring made of a topological su- 
perconductor, there is a MBS on each side of the junc- 
tion and the tunneling current obtains a component that 
is h/e periodic (corresponding to 4tt periodicity in the 
phase difference). This doubling of the flux period with 
respect to the ordinary Josephson effect originates from 
single-electron tunneling mediated by the MBS and is 
dubbed fractional Josephson effect. 

The fo/e-periodic Josephson current is observed as long 
as the fermion number parity of the system is conserved. 
Once the system is in strict thermodynamic equilibrium, 
including relaxation processes which change fermion par- 
ity, the Josephson current reverts to the conventional 
h/2e periodicity. Indeed, the /i/e-periodic Josephson cur- 
rent has equal magnitude but opposite signs for even 
and odd fermion parities, so that it averages to zero in 
the presence of fermion-parity changing processes. Pos- 
sible workarounds that do not require str ict pa rity con- 
servation rely on the ac Josephson effect pZEH Q r finite- 
frequency current noise. 39 Experimental signatures of a 
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fractional Josephson effect in Shapiro step measurements 
have been claimed recently.^ 

Here we show that in mesoscopic rings with a weak 
link, the presence of Majorana fermions can lead to 
an /i/e-periodic Josephson current even in thermody- 
namic equilibrium and in the presence of fermion-parity- 
breaking relaxation processes. This /i/e-periodic contri- 
bution exists in the topological superconducting phase 
and peaks in magnitude near the topological phase tran- 
sition, providing an experimental signature of the phase 
transition. We investigate this signature for a spin- 
less p-wave superconductor wire, the Kitaev chain 
which is a paradigmatic model exhibiting a topological 
phase transition. This model also arises as an effec- 
tive low-energy theory in more realistic situations such 
as the quantum- wire proposals of Refs. [9] and [lOj In 
a ring geometry, the Majorana bound states hybridize 
not only due to the tunneling across the weak link but 
also through the superconducting interior of the ring. 
The latter overlap is exponentially small in the ratio of 
the ring circumference and the superconducting coher- 
ence length governing the spatial extent of the Majo- 
rana bound states. As one approaches the topological 
phase transition, the superconducting coherence length 
diverges and the interior overlap between the Majorana 
bound states becomes significant. This causes a peak of 
the fo/ e- periodic Josephson current near the topological 
phase transition. 

After discussing this effect in clean rings, we extend 
our considerations to disordered rings. We show that 
the signature of the topological phase transition is ro- 
bust and survives under more realistic conditions. This 
issue also leads us to study the influence of disorder in the 
vicinity of the topological phase transition of the Kitaev 
chain which ha d not been discussed previously. Previous 
WQr ]ji6ii 7 | 2Q | 2i | on disorder effects in the Kitaev chain or 
models of quantum wires focused on the regime of large 
chemical potential (measured from the lower band edge), 
\i ^> mA /2 , where A' denotes the effective p-wave order 
parameter of the Kitaev chain in the continuum limit. In 
this regime, the topological region in the phase diagram 
shrinks with increasing disorder.^ In contrast, the topo- 
logical phase transition in the Kitaev chain occurs for 
/jl = and thus in the opposite regime of \i <C mA /2 . Re- 
markably, we find that in this regime disorder increases 
the topological region in the phase diagram. 

This paper is organized as follows. In Section [IT] we 
review the Kitaev model for a one-dimensional spin- 
less p-wave superconductor and its various regimes. We 
also discuss how this model is related to quantum-wire 
based realizations, focusing on the modelling of the mag- 
netic flux through a quantum wire ring in proximity to 
a bulk superconductor. Section |III| is dedicated to the 
flux-periodic Josephson currents in clean rings, focus- 
ing on the flux-periodicity as a signature of the topo- 
logical phase transition. The basic effect is discussed in 



numerical results is given in Sec. |IIID| Sec. |IV| extends 
the considerations to disordered rings. Besides a discus- 
sion of the effects of disorder on the Josephson currents, 
we also study the phase diagram of the disordered wire 
near the topological phase transition. 



II. MODEL 

A. Kitaev model of a one- dimensional spinless 
p-wave superconductor 

Our analysis starts with the Kitaev mo del o f a one- 
dimensional spinless p-wave superconductor^^ 



N 



AT-1 



^ + A TB CjC j+1 + h.c), 



(1) 



which describes a wire of N sites. Electrons on site j are 
annihilated by Cj, hop between neighboring sites with 
hopping amplitude t, and have chemical potential /iTB- 
For all numerical calculations in this paper we choose 
t = 1. The p-wave pairing strength is given by A T b- 
Here, we label both the chemical potential and the pair- 
ing strength by the subscript TB to distinguish these 
parameters of the tight-binding model ([!]) from their 
analogs in the continuum model introduced below. The 
wire can be closed into a ring with a weak link by an 
additional hopping term between sites 1 and TV, 



h.c. 



(2) 



with hopping amplitude t' . 

For an infinite and uniform wire, the Kitaev Hamilto- 
nian exhibits a phase transition when the chemical 
potential /iTB crosses one of the band edges. The sys- 
tem is in a topological (nontopological) superconducting 
phase when the chemical potential is within (outside) 
the interval [— 2t, 2t], i.e., within (outside) the band at 
vanishing pairing Atb = 0- The spectrum exhibits a su- 
perconducting gap on both sides of the topological phase 
transition while the gap closes at the topological critical 
points | /i T B | = 2£. It is thus natural to introduce the 
chemical potential measured from the lower band edge, 
i.e. \i = /i T B + 2£. 

In the vicinity of the band edges (say the lower band 
edge) and thus of the topological phase transition, we can 
make a continuum approximation to the tight-binding 
model ([!]). We will mostly employ the tight-binding 
model in the first part of the manuscript, while we par- 
tially find it more convenient to rely on the continuum ap- 



Sec. IIIB[ analytical considerations on the magnitude of proximation in dealing with effects of disorder in Sec. IV 



the effect are given in Sec. Ill C and a comparison with The continuum model is formulated in terms of the cor- 



3 



(a) 



(b) 




> F A'~ 






i i p 

A' mA' 




(b) 




semiconducting wire 
superconducting ring •* 



FIG. 1. (Color online) Bulk spectrum Eq. Q of Kitaev's 
model for a spinless p-wave superconductor in the regimes (a) 
H > mA' 2 and (b) < \x < mA' 2 . 



responding Bogoliubov-de Gennes HamiltoniaiPE^ 



H 



HA'(x),p} 



(3) 



where A'(x) is the p-wave pairing strength and the curly 
brackets denote the anticommutator. Here, we have in- 
cluded a disorder po tent ial V(x) which we will return to 
in more detail in Sec. IV For V(x) = the bulk spectrum 
of the continuum model is given by 



± 



P 

2m 



|a'|Y 



1/2 



(4) 



which becomes gapless for \i = 0. This indicates the 
above-mentioned topological phase transition between a 
topological phase with ji > and a nontopological phase 
for ji < 0. 

In a semi-infinite wire, the topological phase is charac- 
terized by a Majorana bound state localized near its end 
point. The Majorana bound state has zero energy and a 
wave function that decays exponentially into the wire on 
the scale of the superconducting coherence length £. In 
a finite wire, the Majorana bound states localized at the 
two ends of the wire hybridize and form a conventional 
Dirac fermion whose energy eo scales like the overlap of 
the two Majorana end states which is exponentially small 
in the length L of the wire. The wavefunction of the 
Majorana bound state depends on the parameter regime 
(see, e.g., Ref. 41). This is easily seen by determining the 
allowed wavevectors at zero energy from Eq. Q, which 
yields 



Po 



±zra| A 'l ± V^ra/i - ra 2 |A'| 2 . 



(5) 



(i) ji ^> mA /2 : Deep in the topological phase, the 
bulk excitation spectrum Eq. Q has two minima around 

±Pf = ±y/2mji with a gap « Pf^' (see Fig. Ha). 
According to Eq. ([5|, the Majorana wavefunctions ae- 
cay on the scale £ = 1/raA' and oscillate with a 
much shorter period 1/pf- In a finite wire the hy- 
bridization energy is given by (cf. appendix [A| rp = 
2A'pf\ sm(j9pL)| exp(— L/£), which has accidental de- 
generacies at integer values of pfL/tt. 



FIG. 2. (Color online) Two possible setups for a quantum 
wire with a tunneling junction in proximity to an s-wave su- 
perconductor, (a) The bulk superconductor is interrupted by 
an insulating region underneath the weak link in the wire, (b) 
The bulk superconductor forms a continuous ring and only the 
wire contains a weak link. 



(ii) /i <C mA /2 : Near the topological phase transition 
at fi = 0, the excitation spectrum has only a single mini- 
mum at p = with a gap of order \i (see Fig.J^)). At low 
energies, we can neglect the kinetic energy in Eq. ([3| and 
the spinless p-wave superconductor can be approximately 
described by the Dirac Hamiltonian 



H 



A f pr x . 



(6) 



Eq. gives po ~ iz/i/A 7 , so that the spatial extent 
of the Majorana wavefunction is governed by the co- 
herence length £ = A'//i, which diverges at the topo- 
logical phase transition. In contrast to the previous 
regime, the end-state energy does not exhibit oscillations, 
e ex exp(-L/£). 



B. Magnetic flux 

In the presence of a magnetic flux threading the ring, 
both the tunneling amplitude and the pairing strength 
become complex and acquire a phase. The precise nature 
of these phases depends on the physical realization of 
the Kitaev chain. We illustrate this point by discussing 
two possible setups based on the proposal to realize the 
Kitaev chain in a semiconductor wire proximity coupled 
to an s-wave superconductor as illustrated in Fig. [2J 

(a) The s-wave superconductor is interrupted under- 
neath the weak link in the quantum wire. Current 
can flow around the loop only through the semi- 
conductor weak link. 

(b) The s-wave superconductor forms a closed ring and 
a weak link exists only in the semiconductor wire. 
The current through the weak link of the semicon- 
ductor will in general be only a small perturbation 
of the current flowing through the superconductor. 

We assume that the thickness of the superconducting 
ring is small compared to both its London penetration 
depth and its superconducting coherence length £sc- The 
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supercurrent flowing in the superconductor is given by^ 
2e 

J s = —M 2 (W<p-2eA), (7) 
m* 

where m* and \iJj\ 2 are the effective mass and density of 
the superconducting electrons and ip denotes the phase 
of the s-wave order parameter. The p-wave pairing po- 
tential in the quantum wire inherits its phase cp from the 
s-wave superconductor underneath via the proximity ef- 
fect. (The effective p-wave order parameter may have 
an additional phase shift that depends on geometric de- 
tails such as the direction of the Zeeman field and the 
spin-orbit coupling; however, these contributions lead to 
constant offsets of the phase which are unaffected by the 
magnetic flux.) The vector potential A oriented along 
the wire is related to the Aharonov-Bohm flux <p through 

(j> = j> dxA(x), (8) 

where the integral is taken around the ring of circumfer- 
ence L. 

The phase of the order parameter cp is different for 
the two setups illustrated in Fig. [2] In setup (a), no su- 
percurrent is able to flow since the loop is interrupted, 
J s = 0. If we choose a gauge in which the vector poten- 
tial is uniform around the ring, A(x) = 4>/L, the phase ip 
of the order parameter becomes <p(x) = 47r(0/ '4>q){x / X) 
in terms of the normal- metal flux quantum 0q = h/e. In 
setup (b), the supercurrent around the ring is governed 
by fluxoid quantization cp(x + L) = cp(x) + 27rn, with 
the integer n labeling the fluxoid states. In a gauge in 
which A(x) = 4>/L, this implies that V<p = 27m/L, yield- 
ing a supercurrent of J s = (2e/m*)|^| 2 [27r^n/L — 2eA]. 
Here, [x] denotes the integer closest to x. In thermody- 
namic equilibrium, the system realizes the fluxoid state 
of lowest energy and thus of lowest supercurrent, i.e., 

n=[^/(0o/2)]. 

Within the chosen gauge, in setup (a) the hopping 
amplitude and the pair potential in the tight binding 
Hamiltonian in Eq. M take the form t te i2ncf)/N(f)0 and 
A TB A TB e i47r(</)/ ^^ /Ar) . Alternatively one can per- 
form the gauge transformation Cj — » c ^ e -*O'- 1 / 2 ) 27r 0/ Ar( / ) o 
which eliminates the phase from the pair potential. In 
this new gauge, both the pair potential and the hop- 
ping amplitude t in the interior of the ring are real 
while the hopping amplitude across the weak link ac- 
quires a phase factor, t' —> t'e 1277 ^'^ . Our numerical re- 
sults will be obtained for this representation of the tight- 
binding model. In contrast, in setup (b), we find Atb 
A T Be i27r[0/(</)o/2)](j/Ar) for the pair potential (notice the 
closest integer symbol [.] in the exponent), while t —> 
te i27r0/AT0 o as we n as t t _^ t / e i27r0/AT0 o> As in the previous 

case (a), we can eliminate the phase of the pair poten- 
tial by a gauge transformation. However, this no longer 
eliminates the phase of the hopping matrix element t. 
Instead, one finds t -> te^^^^ ^-^^ ^ and 
if -+ if e iWN){4>/(4>o/2)HN-W/(4>o/2)]} m The fact that we 

can no longer eliminate the magnetic flux from the bulk 



of the wire is a manifestation of the fact that supercur- 
rents in the s-wave superconductor modify the spectrum 
of the quantum wire. 8 

Clearly, the effective Kitaev chain is quite different for 
setups (a) and (b). In the remainder of this paper, we 
will focus on setup (a) where the flux enters only into the 
tunneling Hamiltonian representing the weak link. In this 
setting, the current in the semiconductor wire of interest 
here is experimentally more accessible since there is no 
background current in the bulk s-wave superconductor 
unlike in setup (b). 



III. CLEAN RINGS 
A. Infinite wire 

We first briefly review the Josephson effect of two semi- 
infinite wires connected at their ends through a weak 
link (or equivalently, a ring of infinite circumference), 
as originally considered by KitaevP^ The corresponding 
low-energy excitation spectrum as a function of flux is 
sketched in Fig. [3^. Due to the Majorana end states, 
there are two subgap states whose energies are governed 
by the tunneling amplitude across the weak link. While 
each individual level is periodic in flux with period h/e, 
the overall spectrum is h/2e periodic. As a result, the 
thermodynamic ground state energy of the system - and 
thus the Josephson current in strict thermodynamic equi- 
librium - are h/2e periodic. 

At the same time, the h/e periodicity of the individ- 
ual subgap states is a direct consequence of the Majo- 
rana nature of the endstates. This signature of Majo- 
rana fermions can be brought out in measurements of 
the Josephson current if the fermion parity of the system 
is a good quantum number. The level crossing of the 
two Majorana subgap states in Fig. [3^i is then protected 
by fermion parity conservation. As a result, since there 
is only a single level crossing per superconducting flux 
quantum, the system necessarily goes from the ground 
state to an excited state (or vice versa) when changing 
the flux by h/2e. During this process, the excited state 
is unable to relax to the ground state since this would re- 
quire a change in fermion parity. Thus, the system only 
returns to its initial state after a change in flux of h/e, 
which corresponds to the fractional Josephson effect. 



B. Finite size ring 

For rings with finite circumference, the two Majorana 
bound states localized at the two banks of the weak link 
hybridize not only through the tunnel coupling across 
the weak link but also because of the overlap of their 
wavefunctions in the topological superconductor. In the 
previous subsection, we considered the situation in which 
this interior hybridization is vanishingly small compared 
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FIG. 3. (Color online) (a) Typical Bogoliubov-de Gennes 
spectrum as a function of phase difference across the junction 
between two semi- infinite wires in the topological phase with 
the tunneling amplitude T and the gap A e ff. The two low- 
energy Majorana states represented by the dashed and solid 
lines are related by particle-hole symmetry. The continuum 
of states outside the gap is displayed in gray. The thermody- 
namic ground state has period h/2e. (b) Numerical results 
for the subgap spectrum of a mesoscopic ring with finite cir- 
cumference for A = 1, u = —1.8, t' — 0.01. We set t = 1 
for all numerical calculations in this paper. The parameters 
correspond to £ = 9.5 and the different curves display data 
for ring circumferences L = 95,52,38 all in units of the lat- 
tice spacing. As the circumference of the wire decreases the 
overlap through the topological superconductor in Eq. (12) 
increases. Note that the equilibrium ground state always has 
h/e periodicity in rings of finite circumference. 



to tunneling across the weak link. Conversely, when tun- 
neling across the weak link is negligible compared to the 
interior hybridization, the splitting of the Majoranas due 
to the interior overlap does not depend on flux. Weak 
tunneling across the junction will then cause a small h/e- 
periodic modulation of the split Majorana levels with 
flux. In this situation, even the thermodynamic ground 
state energy becomes h/e periodic, regardless of the pres- 
ence or absence of fermion parity violating processes. In 
fact, of the two /i/e-periodic levels, the negative-energy 
level (which is occupied in equilibrium) corresponds to 
an even-parity ground state while the positive-energy 
level is occupied in the odd-parity first excited state. 
Weak fermion parity violating processes will not destroy 
the /z/e-periodic Josephson current as the two levels no 
longer cross as function of flux. 

The full crossover of the Bogoliubov-de Gennes spec- 
trum as the interior overlap of the Majorana bound states 
increases is illu strated with numerical results in Fig. [3jt 
(see Sec. 



HID 



for details on the numerical calculations). 
They confirm the above picture for the limit of strong 
overlap. But they also show that an /i/e-periodic contri- 
bution to the equilibrium Josephson current exists even 
when the interior splitting is of the order of or smaller 
than the tunnel coupling across the weak link. Indeed, 
the interior overlap essentially pushes one of the two 
states (dashed line) up in energy, while it pushes its 
particle- hole conjugate state (solid line) down. At small 
interior overlaps, this shifts the two level crossings (ini- 
tially at 0o/4 and 30o/4) outwards towards a flux of zero 
and one flux quantum 0q. Note that the level crossings 
remain intact, protected by fermion parity conservation. 
However, once the level crossings reach a flux of zero and 



0o, respectively, the levels merely touch at these points. 
Thus, fermion parity no longer protects the levels from 
splitting, and indeed one state remains at finite and nega- 
tive energies at all values of flux while, symmetrically, its 
particle-hole conjugate state remains at finite and posi- 
tive energies. 

Consider now the Josephson current as function of flux 
in the presence of weak but finite fermion parity violating 
processes. Specifically, we assume that the flux is varied 
by h/e on a time scale which is large compared to the re- 
laxation time of the fermion parity while at the same 
time, the fermion parity violating processes are weak 
compared to the hybridization of the Majorana bound 
states so that the Bogoliubov-de Gennes spectra in Fig. [3] 
are relevant. In this case, the Josephson current is essen- 
tially h/2e periodic deep in the topological phase, where 
L ^> £. However, as the system approaches the topologi- 
cal phase transition, £ grows and hence, the hybridization 
of the Majorana bound states through the interior of the 
ring increases. As a result, the /i/e-periodic contribu- 
tion to the current increases. Conversely, the Majorana 
bound states disappear on the nontopological side of the 
phase transition where the Josephson current thus re- 
verts to h/2e periodicity. As a result, we expect a peak in 
the fo/e-periodic Josephson current near the topological 
phase transition, whose measurement would constitute 
a clear signature of the topological phase transition and 
the formation of Majorana fermions. 

This expectation is confirmed by the numerical results 
shown in Fig. [4^i, where the corresponding Fourier coeffi- 
cient A h j e = (2e/h) f^ e dcj) I ((/)) sin(27r ecj)/h) of the equi- 
librium Josephson current I{(j)) is plotted as a function of 
chemical potential. A h / e ([i) exhibits a peak in the topo- 
logical phase (/i > 0), which moves closer to the topolog- 
ical phase transition at \i = as the ring circumference 
increases (see Fig. PI 



Deep in the topological phase the Majorana bound 
states are localized at the weak link. Approaching the 
phase transition at \i = 0, the MBS delocalize. On the 
one hand, this causes an increase in the overlap of the 
MBS in the interior of the topological superconductor. 
As discussed above, this leads to an increase of the h/e- 
periodic Josephson current. On the other hand, however, 
the probability density of the Majorana bound state near 
the weak link decreases, causing a suppression of the hy- 
bridization of the Majorana bound states across the weak 
link and hence of the /i/e-periodic Josephson current. 
Thus, the peak occurs for the value of \i where the in- 
terior overlap splitting is equal to the tunnel coupling. 
Since the interior overlap is exponentially small in L/£ 
while the hybridization across the weak link is roughly 
independent of the ring's circumference L, the peak po- 
sition shifts towards the phase transition point at \i = 
with increasing L (cf. Fig. [4^). Since at the same time 
the /z/e-periodic Josephson current becomes suppressed 
when the systems is approaching the phase transition, 
the peak is more pronounced in shorter rings. 
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C. Low-energy Hamiltonian 



(a) 



A more quantitative description can be developed by 
restricting the Hamiltonian to the low-energy subspace 
spanned by the two Majorana bound states. The projec- 
tion of the tunneling Hamiltonian across the weak link 
onto this subspace yields 



H T = -rcos(27T0/0o)(dL d M - 1/2), 



(9) 



where c?m is the Dirac fermion constructed from the two 
Majorana bound states. The parameter T measures the 
tunnel coupling of the Majorana b ound states across the 
weak link and is given by (cf. Eq. (A4) in appendix |A|) 



t>(4t - /i)A T B 
t(t + Atb) 2 " 



(10) 



Here the factor of \i accounts for the probability den- 
sity of the Majorana wavefunction at the junction, which 
vanishes at the phase transition. 

The overlap of the Majorana end-states in the inte- 
rior of the wire leads to an additional coupling (cf. ap- 
pendix [a| 



overlap — ^0 



(d^d M ~ 1/2) , 



where 



e = 2/iexp(-L/£) 



(11) 



measures the strength of the overlap. 

Combining these two contributions for a mesoscopic 
ring near the topological phase transition (/i <C mA /2 ), 
the effective low-energy Hamiltonian reads as 



rcos hd 



(4d M - 1/2) . (12) 



The Bogoliubov-de Gennes spectrum of this Hamiltonian 
reproduces the numerically calculated subgap spectra de- 
picted in Fig. [3)3. 

In principle, both the negative energy continuum states 
as well as the negative energy subgap state contribute 
to the equilibrium Josephson current. If we denote the 
sum over all negative excitation energies by Eq ((/)), we 
can write the equilibrium Josephson current as I{(j)) = 
—d^Eo^cj)). However, it is natural to expect and will be 
coroborated by our numerical results that the Joseph- 
son current is dominated by the contribution of the sub- 
gap state 7(0) ~ d^eo — T cos(27T0/(/>o)|/2. Thus, it 
is straight-forward to compute the /i/e-periodic Fourier 
component of the Josephson current, 



e r 



1 - ^ -f arcsin(f ] 



2/1' 



, e <r 
60 >r 



. (13) 



In the next section, we compare this analytical result 
with numerics and find nice agreement. 
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FIG. 4. (Color online) (a) Numerical results for the h/e- 
periodic Fourier component of the Josephson current, A h / e , 
as a function of chemical potential in a ring with Atb = 1 
and L — 200 (blue dots) together with analytical expression 
Eqs. (13) (red solid line). Inset: numerical results for eo (blue 
squares) and F (red circles) together with the correspond- 
ing analytical expressions (gray dashed curves) Eqs. ( |11| ) 
and (10). (b) /i/e-periodic Fourier component (solid), h/2e- 
periodic Fourier component (dashed) , and the maximum tun- 



neling current of the MBS, eT/h. 
circumferences L. 



(c) A h / e for different ring 
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D. Numerical Results 



To obtain numerical results for the Josephson current, 
we solve the Hamiltonian defined in Eqs. ([I]) and Q by 
exact diagonalization. Fig. [4^i compares the amplitude 
of the /i/e-periodic component as a function of chemical 
potential with the analytical result in Eq. (13). The nu- 



merical results agree well with the behavior predicted by 
the low-energy model, except for small deviations in the 
immediate vicinity of the phase transition at \i = 0. In 
the inset of Fig. [4^i we compare the analytical and nu- 
merical results for the quantities T and eo appearing in 
the low-energy Hamiltonian. While the model correctly 
captures eo in the regime of interest, there are devia- 
tions of T near \i = 0. These discrepancies are readily 
understood as a consequence of the finite circumference 
of the ring. Although the coherence length diverges at 
the phase transition, the Majorana bound states can de- 
localize at most throughout the entire length of the ring 
there remains a finite probability density of the Majorana 
bound state wavefunction at the weak link. 

Figure [4b shows that the left flank of the peak of A h j e 
and eT / //"deviate slightly. This deviation is a measure of 
the size of the bulk contribution to the /i/e-periodic cur- 
rent. The latter can thus be seen to be small, justifying 
our focus on the low-energy Hamiltonian (12) describing 
the Majorana bound states only. In the same figure, the 
h/2e component is plotted, showing that the /i/e-periodic 
Josephson current exceeds the h/2e component. This is 
a consequence of the tunneling regime that favors single- 
electron tunneling over the tunneling of Cooper pairs. 

In Fig. 4c, we show how the position of the peak in 
Ah/e depends on the circumference of the ring. We find 
that the value of \i where the peak occurs scales as 1/L. 
This result can be understood as follows. T is essentially 
independent of the length of the ring, while eo scales as 
~ exp(— L/£). As we have seen above the peak occurs 
at eo = r. For given t, Atb> and t' ', V is fixed and the 
peak occurs at a constant value of the ratio L/%. Since 
£ ~ l//i, the value of /i where the peak occurs scales as 
1/L. Also note that the above-mentioned tail of the peak 
at \i < originating from finite-size corrections is more 
pronounced in shorter rings. 



(a) A h/e (t'c/fl) 
0.25 






(b) A h/e {Vm 
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(c) A h/e (Ve/fi) 
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IV. EFFECTS OF DISORDER 

A. /i/e-periodic Josephson current in disordered 
rings 

In this section we investigate the fate of the peak in the 
equilibrium /i/e-periodic Josephson current in the pres- 
ence of disorder. Our main results are: 

(i) The typical peak height is not affected by disorder 
as long as the mean free path is longer than the cir- 
cumference of the ring. Thus the signature persists 
in the presence of moderate disorder. 



FIG. 5. (Color online) (a) Numerical results for the h/e- 
periodic Fourier component of the Josephson current, A h / e , 
as function of chemical potential for a clean ring (black solid 
line) and disordered rings with four disorder configurations 
(green dashed lines) corresponding to / = 5. (b) A h / e for a 
clean wire (black solid line) together with the histogram of 
the peak position in the presence of disorder for / = 75 as a 
color code (green (gray) area), (c) Same as (b) with / = 5. 
For all plots we chose L = 20 and Atb = 1. 
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0.4 0.5 r A' 0.25 0.5 0.75 r/t* 



FIG. 6. (Color online) Histogram of F for the same parame- 
ters as in Figs. [5|) and c at /x = 0.4. The dashed line denotes 
the value of F for the clean ring. 
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(ii) For stronger disorder the peak height decreases and 
the peak position is shifted to lower chemical po- 
tentials. 

To study the effect of disorder we add a random on- 
site potential V^cJq to the tight-binding Hamilto- 
nian and Q, where the Vi are taken from a uni- 
form distribution over the interval [-W, +W]. The mean 
free path is then related to the disorder strength as 
I oc 1/VF 2 P^ To obtain numerical results we compute 
the spectrum by exact diagonalization. Disorder affects 
the /z/e-periodic Josephson current by introducing fluc- 
tuations in the quantities eo and F. While F is mainly 
affected by local fluctuations of the probability density 
of the Majorana wavefunction at the junction, eo fluctu- 
ates due to the disorder potential in the entire ring. The 
interior overlap in disordered wires has been investigated 
previously for the continuum model ([3| in regime (i), i.e., 
[i ^> m A' 2 pH where disorder leads to an increase of eo 
and subsequently to a disorder-induced phase transition 
to the nontopological phase. 

Fig. [5^1 shows numerical results for the /i/e-periodic 
Josephson current for a few disorder configurations. The 
peaks in the presence of disorder (green dashed curves) 
are of comparable height as the peak in the clean ring 
(black solid curve). The peak shifts as a function of 
chemical potential which indicates fluctuations of the co- 
herence length due to disorder. 

Surprisingly, the peak shifts to lower chemical poten- 
tials, corresponding to a decrease in eo with disorder in 
stark contrast to the known case of large \i. This implies 
that the topological phase is stabilized by disorder if the 
system is close to the phase transition. To investigate 
this further we plot the height and position of the peak 
maxima of many disorder configurations as a color code 
histogram for I > L in Fig. [5]o and I < L in Fig. [5J3. In- 
deed the average peak height is comparable to the one in 
the clean case for I > L. When I < L the average peak 
height starts to decrease. The histogram in Fig. ^ con- 
firms that the peak is shifted to lower chemical potentials 
on average. 

To understand this behavior we analyze the probabil- 
ity distributions of eo and F over the disorder ensemble. 
In Fig. [6] we show numerical results for the histogram of 
T corresponding to the two ensembles in Figs. [5b and c at 



FIG. 7. (Color online) Slope of (ln(e (L)/2A e fr)) (see inset) 
vs. disorder strength l/l for regimes (i) with \i — 300 mA /2 
(blue dots) and (ii) with \i — 3 x 10 -3 mA /2 (red cros ses) 
together with theoretical prediction according to Eqs. (14) 
(dashed line) and (15) with A = 1/2 (solid line). Inset: nu- 
merical data (red crosses) and linear fit (solid) of the average 
of In eo as a function of L for different disorder strengths. 



H = 0.4. For weak disorder the distribution is symmet- 
ric with a mean near the zero-disorder tunnel coupling. 
For larger disorder when I < L the distribution becomes 
wider and asymmetric and the average decreases. 

In order to determine the probability distribution for eo 
we now turn to the continuum Hamiltonian ([3| for a wire 
of length L without tunnel junction. To model short- 
range correlated disorder in the continuum model, we 
include a disorder potential with zero average (V(x)) = 
and correlation function (V(x)V(x')} — 7$ (a; — x'). For 
this model we employ a numerical method based on a 
scattering matrix approach! 20 * 44 * 45 ! From the scattering 
matrix S we obtain the lowest energy eigenstate eo by 
finding the roots of det(l — S(e)). In this model, the 
probability distribution of the hybridization energy eo has 
been shown to be log-normal in Ref. |2lJ Specifically, it 
was shown that the log-normal distribution is governed 

by 



(ln(e /2A eff )) = -L 



£ 21 



varln (e /2A eff ) = 



L 
21' 



(14) 



for regime (i). The distribution function reflects the 
disorder-induced phase transition to the nontopological 
state at £ = 21. 

The numerical results are presented in Fig. [7[ In the 
inset, we show that the mean of ln(eo/2A e ff) is indeed lin- 
ear in L with the slope depending on disorder strength. 
This slope is plotted as a function of inverse mean free 
path in Fig. [7| The data for fi = 300mA /2 (blue dots) 



agrees well with the prediction Eq. (14) with the defini- 
tions I = Vjp/j and £ = 1/raA'. 

The same plot also shows data corresponding to regime 
(ii), i.e., ijl <C mA /2 , marked by red crosses. Here, we 
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have I = A /2 /7 and £ = A'//i. Clearly, the behavior 
is qualitatively different from regime (i), since disorder 
decreases eo rather than increasing it. This is consistent 
with the shift of the peak of A h j e to lower \i. 

In order to gain analytical insight we now derive the 
probability distribution of eo in regime (ii) extending the 
results of Ref. '2TJ The relevant momenta at low energies 
in this regime are near p = (cf. Fig.[l]3). Linearizing the 
dispersion around this point yields the Dirac Hamiltonian 
Eq. ^ , where the disorder potential enters as a random 
mass term. Since the disorder potential is short-range 
correlated it couples high- and low-momentum degrees 
of freedom in the original Hamiltonian. Thus a proper 
linearization of the Hamiltonian requires one to project 
out the high-momentum states, which renormalizes the 
gap. 

For a strictly linear model with a random mass term, 
the overlap eo has a log-normal distribution,^ 



(ln(e /2A eff )) 
var ln(e /2A eff ) 



L 

T 



i 



Thus for the Dirac Hamiltonian the mean of In (eo/2A e ff ) 
does not depend on disorder. A systematic lineariza- 
tion of the disordered spinless p-wave superconductor in 
the vicinity of the topological phase transition effectively 
renormalizes the chemical potential /i and hence the co- 
herence length £ = A'//i. 

We start by defining the projection operators P = 
^\p\<pi I^p) onto ^ e l° w momentum subspace and 
Q = 1 — P, where {\^ p )} p is a complete set of momentum 
eigenstates of the clean Hamiltonian. The relevant mo- 
mentum scale for this projection is given by pi = mA', 
since for p <C pi, the term p 2 /2m constitutes the lowest 
energy scale of the Kitaev Hamiltonian. Furthermore, we 
assume that the disorder potential does not affect high 
momenta p\ ^ 1/1. We can now project the clean Kitaev 
Hamiltonian H to the low- and high-energy subspaces, 



PHP~P[(-v + V(x))t z 
QHQ ~ Q (p 2 /2m) r z Q. 



■ A 'pr x \ P, 



(16) 
(17) 



Both subspaces are exclusively mixed by the disorder po- 
tential PHQ = PV (x)r z Q . To second order in V, the 
correction to the low-energy Hamiltonian is then given 

by 



8H(p) ~ | PHQ (e p - QHQ)' 1 QHP V P 

1 



E v f 



pp' 



\p'\>pi 



e p — p' 2 /2mr 2 



V n >r 



(18) 



Here we used the short notation V vv > = (ip p | V(x) \ ip p r 
Averaging over disorder, we obtain 



(SH(p)) 



E 

\p'\>pi 



2m r ) 
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(19) 



Thus the renormalization produces a contribution to 
the low-energy Hamiltonian which has the same struc- 
ture as the chemical potential term. Hence we find a 
renormalized chemical potential \jl = \i + \j/A' with 
a numerical factor A > that cannot be determined 
from this argument. Thus disorder enters the final re- 
sult through the renormalized coherence length 



A 



1 1 



(20) 



The data in Fig. [7] confirm Eqs. (15) and (20) and de- 
termine the unknown numerical prefactor to be A = 1/2. 
Thus for \i <C mA /2 , eo has a log- normal distribution 
with mean and variance given by 



(ln(e /2A eff )) = -L 
varln(e /2A eff ) = y. 



1 1 

£ + 21 



(21) 



(15) 

This result is very similar to Eq. fll4fo where, however, the 



disorder correction to the decay length enters with oppo- 
site sign. This underlines the contrast between the two 
regimes, i.e., that disorder drives the system further into 
the topological phase when it is close to the phase tran- 
sition, but away from it for larger chemical potentials. 
Specifically a spinless p-wave superconducting wire with 
negative chemical potential may exhibit edge states with 
an energy exponentially small in L as long as disorder is 
strong enough. 

Combining the disorder-induced fluctuations of T and 
eo we can understand the suppression of the peak in 
the /i/e-periodic Josephson current in Fig. [5^ for I < L. 
While eo is decreased on average for a given \i with in- 
creasing disorder, T does not increase at the same time 
and thus the average peak height decreases. However 
the fluctuations of T and eo become larger as disorder 
increases such that for single disorder configurations sig- 
nificant peaks are still possible even if the average peak 
height decreases. 



B. Phase diagram of a disordered wire 

Motivated by the contrasting probability distributions 
of eo in the regimes of large and small \i we numerically 
calculate the phase diagram of the continuum model (|3| 
as a function of /i and 7, particularly paying attention 
to the region near the topological phase transition of the 
clean model. By means of the scattering matrix approach 
also used in the last section we compute the determinant 
of the reflection matrix of a wire of length L at e = 
which approaches the values +1 and — 1 as L 00 i n th e 
nontopological and topological phase, respectively! 32 * 46 * 

The resulting phase diagram is plotted in Fig. [8] From 
Eq. (14) we infer that the topological phase transition 
occurs for £ = 21 in the regime \i ^> mA /2 . Using the 



definitions of / and £ in this regime, we obtain the phase 
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u <D 



H (mA' 2 ) 



FIG. 8. (Color online) Phase diagram of the continuum model 
|3| as function of disorder strength 7 and chemical potential 
\i in the regime \i <C mA /2 . The data has been averaged over 
100 disorder configurations. For \i < disorder gives rise 
to a trivial-to-topological phase transition with a reentrant 
nontopological phase for stronger disorder. The dashed line 
denotes the pha se transition line 7c*^(/i) valid for small 
given in Eq. |22| . Inset: Phase diagram for a larger range of \i 
and 7. The solid line represents the predicted phase boundary 
7c^(/i) for large \i. (The analytical phase boundary is only 
accurate at large [i up to sublinear corrections.) 



boundary Jc\/jl) = 4A'/i. This is compared with the 
numerical results in the inset of Fig. [8] The numeri- 
cally calculated phase boundary 7c UIT1 (aO has only sub- 
linear deviations from the predicted line, so that the ratio 
7j mm (/i)/2/iA / = £/Z approaches the value 2 for fi — >> 00 
as expected. 

However, near \i = the behavior is qualitatively dif- 
ferent. Here, disorder can induce a topological phase 
for ji < as well as a reentrant nontopological phase 
at larger disorder. From Eq. (21) we find the condition 
£ = — 21 for the phase boundary. This corresponds to 



decreases again. 



V. CONCLUSION 

Even for conventional superconducting phases, the flux 
periodic currents have been widely studied for meso- 
scopic rings we studied the Josephson currents 

across a weak link in a mesoscopic ring in a topological su- 
perconducting phase. As a paradigmatic model system, 
we studied Kitaev's model of a one-dimensional spin- 
less p-wave superconductor, focusing on the parameter 
regime near the topological phase transition. We found 
that in mesoscopic rings, there is an /i/e-periodic contri- 
bution to the tunneling current even if electron number 
parity is not conserved. This /i/e-periodic contribution 
emerges due to the hybridization of the Major ana bound 
states localized on the two sides of the weak link through 
the interior of the ring and exhibits a pronounced peak 
just on the topological side of the topological phase tran- 
sition. This peak provides an interesting signature for 
the existence of a topological phase transition and the 
formation of Majorana fermions at the junction. 

We found that this effect remains robust in the pres- 
ence of disorder in the wire. In fact, near the topological 
phase transition disorder can even stabilize the topolog- 
ical phase. When tuning, say, the chemical potential of 
the system to the nontopological side of the phase tran- 
sition, there is a disorder-induced topological phase for 
moderate amounts of disorder, with a reentrant nontopo- 
logical phase at even stronger disorder. This is in stark 
contrast to the behavior of the system far in the topolog- 
ical phase where disorder weakens and eventually desta- 
bilizes the topological phase. 
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7 f ) (M) = -2A / /i, 



(22) 



which we find to agree well with the numerical results 
for the phase diagram (see dashed line in Fig. [8|. Thus 
the phase diagram confirms that weak disorder leads to 
an enhancement of the chemical potential range of the 
topological phase, while for stronger disorder the range 
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Appendix A: Effective energy splitting and Josephson coupling of Majorana end-states 



We present here the analytical estimation of the Majorana energy splitting, eo, and the effective Josephson coupling, 



T, in the Hamiltonian (12). For sake of simplicity we will first present our calculation for the effective Josephson 
coupling T, obtained by working in the tight-binding model. We will then describe the calculation of eo through 
an alternative approach working directly in the continuum limit both in both regimes \i <C mA' 2 and fi ^> mA' 2 
discussed in the main text. 

Effective Josephson coupling. In order to compute the effective Josephson coupling we neglect the effect of 
the overlap of the two Majorana wavefunctions in the topological part of the ring. We therefore consider the limit 
of a junction between two half-infinite topological sectors of the wire, the right one on sites j G [l,oo) and the left 
one on j G (oo, —1). They are both described by the Hamiltonian with paring strength AxBe^ a with a = L,R 
for the left and right sector, respectively. The hopping between the two sectors from Eq. ([2| now simply reads 
H T = -t'(cl lCl +c t 1 c_i)— c.f. Fig.^ 

We include the effect of the tunneling between the two wires perturbatively. The low energy excitations of this 
model for t' = are represented by left and right zero energy Majorana operators^ 

oo 

bR = Aj2(*i-x j -)l ( B R l (Al) 

— oo 

b L = A aC'hS (A2) 

j=-i 

(A3) 

where x± (— /xtb ± v//4 B -4t 2 +4A* B )/(2i + 2A TB ), A = (£ | x\ — x J _\ 2 ) 1 ' 2 is a normalization constant, and 

the operators Ja(b)j are defined via cj = (7^e _ *^ R / 2 + i^j^e l ^ R / 2 ) /2. The projection of onto the subspace 
spanned by the operators 6r, b^ leads to the effective coupling between the Majorana states. The Hamiltonian can 
be rewritten in terms of ordinary fermion operators dM = (&r + z&z,)/2 to take the form presented in Eq. Q with 

T = t' A 2 \x+ — X-\ 2 . For simplicity, we consider fi < 2t (l — ^/l — A^ B /t 2 ^ which corresponds to the condition 

H < mA ' 2 in the continuum model ^ and ensures that x± are real (see discussion of the continuum model in 



Sec. II A). Explicitly one obtains 



^(4t-^)A TB 

tit + Atb) 2 ( ' 

In the continuum limit, when \i <C Atb <^ t the expressions simplifies to 

T « 4t'iiA TB /t 2 . (A5) 

Majorana energy splitting. In order to compute the energy splitting eo we employ an alternative method 
working directly in the continuum limit. Similar to before, we neglect here the Majoranas interaction through the 
Josephson junction. The problem is then completely equivalent to calculating the energy splitting of two Majorana 
at the end of a wire of length L. We start considering the regime (ii) described by the Hamiltonian in Eq. (|6|, 
H = h(x)t z + A'(x)pr x , with the specific choice of parameters (c.f. Fig. [9]) 

fj,(x) = -V [O(-x) + e(x - L)] + fj,&(x)e{L - x) , 

A\x) = A' , (A6) 

where Vb > and \i > guarantee that the wire is in a topological phase in [0, L] and in a nontopological phase 
otherwise. 

We determine the Majorana wavefunctions by a perturbative approach. We first find the wavefunctions of the 
Majorana fermion localized at one of the interfaces between the topological and the insulating region, thus fully 
ignoring the existence of the other Majorana state. We label the corresponding states at the interfaces at x = 
and x = L as \L) and \R) respectively. We then project the Hamiltonian onto the Majorana subspace, to obtain 
their effective interaction, i^overiap = ^2 a p=L R( a \H\P) \ a ) Eventually we will be interested in the limit Vb — » oo 
corresponding to a high insulating barrier outside the wire. 

The Majorana state at x = is the zero-energy eigenstate of the Hamiltonian H^, defined again by the same 
Hamiltonian as in Eq. |6]), but with the choice of parameters 

fi(x) = -VqG(-x) + fiS(x) , 

A'(x) = A' , (A7) 
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FIG. 9. (a) Sketch of the Josephson junction between two topological segments of the wire and spatial dependence of the 
gate potential in the corresponding continuum realization, (b) Spatial dependence of the chemical potential, in the 

Hamiltonian of the finite length wire and for the approximate Hamiltonians, Hl and Hr, used in the perturbative calculations 



as depicted in Fig. [9j Solving this equation separately for x > and x < leads to zero energy states with imaginary 
momenta. Namely we can write 



\L) = v_, L e aoX e(-x) + v + , L e- ax G(x) 



(A8) 



where a = fi/\A'\ is the inverse coherence length in the wire, and ao = Vo/\A'\. The twodimensional vectors v±,l 
are determined by the continuity of the wavefunction and its derivative at the interface and by the wavefunction 
normalization. They are given by: 



V+,L = V_ L 

In the limit Vq — >• oo we are interested in, they reduce to 



a + ao 



1 



(A9) 



(A10) 



In a completely analogous way we can calculate the zero-energy eigenstate of the Hamiltonian Hf> defined once more 
by the Hamiltonian in Eq. (J6|, now with (c.f. Fig. [9]) 

H(x) = -V O(x - L) + fiQ(L - x) , (All) 
A'(x) = A' . (A12) 

In this case the zero-energy eigenstate reads 



\R) = 



a + ao 



-OLo{x-L)r\ 



G(x - L) + e' 



0(L - x) 



(A13) 



In the limit Vq — » oo the prefactor reduces to y/a. Note that the particle-hole superposition has different phases in 
\L) and \R). 

We can then project the full Hamiltonian onto the low-energy subspace spanned by the two Majorana states. 
In doing so, we conveniently rewrite H = Hl + Vr = Hr + Vl, where Vr = — (Vb + /jl)Q(x — L)r z and Vl = 
— (Vb + /i)6(— x)r 2 , and the spatial dependence of the various terms is presented in Fig. [9] We can then compute, e.g., 
(L\H\L) = (L\Hl\L) + (L\Vr\L) = (L\Vr\L), and all the other matrix elements in a similar way, to get 



_ f(L\V R \L) (L\V L \R)\ 

^overlap- I {R \ Vl \ L ) (R\ Vr \R) ~ 











(A14) 



This leads to the energy splitting eo = 2/ie aL , as reported in Section III C 



In complete analogy, one can perform the calculation of the energy splitting in the regime \i ^> mA' 2 . In this case 
the appropriate Hamiltonian, H, is that in Eq. ([3|, with the same choice of parameters as in Eq. (A6). Again we 
start considering the left interface, looking for zero-energy solutions of the Hl, i.e., the Hamiltonian in Eq. Q with 
the choice of parameters as in Eq. (A7). We introduce ao = Vq/ {mA' 2 ) > and a = /i/(mA /2 ) > 0. We can write 
the solutions as 



\L) = vW^^ ^ [(1 + ri)e iK -*Q{-x) + (e ifc +* + V e ik ~ x ) Q(x)] , 
\R) = vWj VT~^T (]) [(e- ik +( x -V + T]e- ik -^ x - L ^ 9(L - x) + (1 + r])e- iK - x O(a 



l + 2a 



L) 



(A15) 
(A16) 
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where k>- = — im\A f \^ 



k± = m\A , \(±y / 2a + i) and r] — (k,- — fc+)/(fc_ — «_). In all the expressions we are 



neglecting O(y/a/y/ao, l/^/oo), consistent with the condidtion ao ^> a ^> 1, reflecting the limit - 
under consideration. The matrix elements of the effective Hamiltonian are, in this case, (L\H\L) 
(L\H\R) = (R\H\L) = -c , with 



oo and the regime 
: (R\H\R) = and 



eo = A' ^/2mfi sin(^2ra/iL) — ^/^/^(^ — 1) cos(v // 2m/i)L + 0(y/a/afy 



(A17) 



